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I. INTRODUCTION 



Key electroweak observables of the standard model (SM), such as the mass and width 
of W bosons, will be measured with high precision in production of W"^ and bosons in 
hadron collisions at the Fermilab (Tevatron) and Large Hadron Collider at CERN (LHC). 
The uncertainty in the W boson mass Mw will be reduced to 30-40 MeV per experiment in 
the Tevatron Run-2 and to about 15 MeV at the LHC At this level of accuracy, the 

theory framework must incorporate effects were dismissed in less precise analyses. 

One such feature is the dependence on masses of quarks, which is frequently neglected 
in hard-scattering reactions. While the massless approximation is clearly appropriate in 
channels involving only light quarks (m, d, and perhaps s), it may be less adequate when 
heavy quarks (c and b) are involved, particularly when the cross section depends on a small 
momentum scale close to the mass mg of the heavy quark. The distributions of heavy 
electroweak bosons over their transverse momenta qr may be sensitive to the masses of c 
and b quarks, given that the transverse momenta of most bosons (of order, or less than a 
few GeV) are comparable to uiq. The impact of the quark masses on the qr distributions of 
bosons with invariant masses Q > Mw is suppressed by properties of soft parton radiation, 
as demonstrated below. Nonetheless, the quark mass effects are relevant in high-precision 
studies, such as the measurement of M^, or in reactions dominated by scattering of heavy 
quarks (particularly bottom quarks), such as Higgs boson production in bb annihilation. 

In the present paper, we quantify the effects of the heavy-quark masses on g-r distributions 
in Drell-Yan production of W^, Z^, and Higgs bosons. We evaluate the associated impact 
on the W boson mass measurement and demonstrate the importance of the quark-mass 
corrections in the processes initiated by bottom quarks, like bb — > Higgs. This latter channel 
can have a large cross section in supersymmetric extensions of the standard model 0,0,01. 

The description of the heavy-quark scattering is complex because of the presence of the 
heavy-quark mass scale mg, in addition to the boson's transverse momentum and its vir- 
tualit y Q . Two popular factorization schemes (fixed-flavor number (FFN) scheme 0, Q, 0, 
and zero-mass variable flavor number (ZM-VFN) scheme) were applied recently to 
calculate various observables in the 6-quark channels |12 . How- 

ever, neither of the two schemes is entirely consistent when describing logarithmic corrections 
in processes like bb ^ H in the small-gr region. When is much smaller than Q, the calcula- 
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tion of the qx distribution in perturbative quantum chromodynamics (PQCD) must evaluate 
an all-order sum of large logarithms ln(g|,/Q^), arising as a consequence of the recoil of the 
electroweak bosons against unobserved hadrons of relatively low energy or transverse momen- 
tum (soft and collinear hadrons) . The resummation of Drell-Yan qr di stributions for massless 
initial-state quarks has been developed in a variety of forms Q Q H Q, Q- In 

the present paper we use the Collins-Soper-Sterman formalism (CSS) (2J|, which excellently 
describes the available gr data from fixed-target Drell-Yan experiments and boson pro- 
duction at the Tevatron. When applied to heavy-quark scattering, small-gr resummation 
indispensably involves logarithms depending on both the transverse momentum, ln(g|./(5^), 
and heavy-quark mass, ln(/i|./mg). To correctly treat both types of large contributions, 
Ref. 3l 

extended the CSS formalism to a general-mass variable flavor number (GM-VFN) 
scheme 2^- The extended resummation formalism was then applied to describe heavy- flavor 
production in deep inelastic scattering (DIS). We review the qx resummation formalism for 
heavy quarks in Section ITTl 

The present study was motivated by the results of Ref. where the correct treatment 
of the heavy-quark masses was found to substantially change the differential cross sections 
for heavy quark production in DIS at small qx ~ mg. In Section lTTT| we extend that study to 
examine the heavy-flavor effects on g-r distributions of W^, and Higgs bosons produced 
at the Tevatron and LHC. To get a flrst idea about the magnitude of the mass effects, we 
consider fractional contributions of reaction channels with initial-state c and h quarks to the 
entire production rate (see Table P). At the Tevatron, the heavy-quark channels contribute 
only 8% (3%) of the inclusive cross section in (Z°) boson production, and consequently 
the quark masses can be usually neglected. At the LHC, the heavy-quark channels add up 
to 22% in and 31% in W~ boson production. Modiflcations caused by the heavy-quark 
masses at the LHC are comparable to the other uncertainties and must be considered in 
precision measurements of Mw- We then turn to Higgs boson production in the bb H 
channel. We show that the nonzero mass of the bottom quark substantially modifies the 
small-gT Higgs distribution. Furthermore, we note that other processes with initial-state 
heavy quarks, like bg Z^b, H^t, H^i, H^b, etc. will be affected by the heavy-quark 
masses in a similar way. We summarize the results of this study in Section Hvl 
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II. TRANSVERSE MOMENTUM RESUMMATION FOR MASSIVE QUARKS 



A. General form of the resummed form factor 

In this section, we define the essential elements for the transverse momentum resum- 
mation in the presence of massive quarks. The resummed differential cross section for the 
inclusive production of heavy electroweak bosons in scattering of initial-state hadrons A and 
B takes the form 
da r bdb 



dQ'^dydqj 



Q 

bdb ~ 

/ — MqTb)W{b,Q,XA,XB,{mQ}) + Y{qT,Q,y,{mQ}), (1) 

dT Jo 

where y = (1/2) In [{E + pz) /{E — p^)] is the rapidity of the vector boson, xa,b = 
are the Born-level partonic momentum fractions, S is the square of the center-of-mass energy 
of the collider, and Joiqrb) is the Bessel function. The integral is the Fourier-Bessel transform 
of a form factor Q, a;^, a;^, {m^}) in impact parameter (6) space, which contains the 
all-order sum of the logarithms ln'"(g|,/(5^). The 6-space form factor is given by 

W{b,Q,XA,XBd^Q}) = ^ 5^ ag) e-^^^'Q'^'"^-^" Vj/a{xaA {^q}) ^/bIx,,, 6, W), (2) 

where the summation is performed over the relevant parton fiavors j and k. Here, aj^'' is the 

product of the Born-level prefactors, e"'^^''''^'^'"^^^^ is an exponential of the Sudakov integral 

/■^' du^ r /OS 1 

Sib,Q,{mQ}) = / ^ Aia,ifi),{mQ})ln(^)+Biasm)AmQ}) , (3) 

and Vj/A{x,b, {itiq}) are the 6-dependent parton distributions. The constant factor 6o = 
2e~'^'^ ~ 1.123 appears in several places when a momentum scale is constructed from the 
impact parameter, as in the lower limit b^/b"^ of integration in Eq. The Sudakov expo- 
nential and 6-dependent parton densities resum contributions from soft and collinear multi- 
parton radiation, respectively. In the perturbative region (6^ <^ Ag^^), the distributions 
'Pj/A{x,b, {mQ}) factorize as 

— Cj/a{x/i,b, {mgj./ii.) fa/A{^,f^F) 
a=g,u,a,... ^ 

= {Cj/a ® fa/ a) (x, b, {tUq}, ^p) (4) 

into a sum of convolutions of the Wilson coefficient functions Cj/a{x,b, {fnQ},fiF) and k^- 
integrated parton distributions fa/A^i, I^f)- Note that the initial state A can be a hadron as 
well as a parton, in which case we denote the partonic initial state as A' . 



If quark masses {ruQ} are neglected in the coefficient function Cj/a{x,b, {itiq}, fip), it 
depends on b and fip only through the logarithm In^fipb/bo): 

lim Cj/aix,b,{mQ}, jip) = Cj/aix,\n^^). (5) 

{bmQ}^{0} Do 

For this reason, we set the factorization scale fip equal to bo/b to prevent the large logarithms 
ln(/i^6/6o) from appearing in Cj/a in the limit 6 — > 0. 

'^{Qt, Q, y, {^q}) in Eq. ^ is the regular part, defined as the difference of the fixed-order 
cross section and the expansion of the Fourier-Bessel integral to the same order of it 
dominates at gr ~ <5 and is small at — > 0. The regular piece Y and the dominant contri- 
butions in W can be calculated in PQCD, if Q is sufficiently large. A small nonperturbative 
component in W contributing at Qt less than a few GeV can be approximated within one 
of the available models. Our specific choice of the nonperturbative model is described in 
Section mn 



B. Extension to the case of massive quarks 

Eqs. Ill) and Q indicate that W and Y may explicitly depend on the masses of heavy 
quarks {rriQ}, with Q = c, b, and t. Masses of the light u, d, and s quarks {mu^d,s ^ 
Aqcd) are neglected by definition in hard-scattering contributions and implicitly retained 
in the nonperturbative functions, i.e., parton densities fa/A{x,fiF) and power corrections 
to W at large impact parameters b > I GeV~^. Indeed, the nonperturbative dynamics 
described by fa/A{x,fip) and power corrections does depend on the light-quark masses, but 
this dependence is not separated explicitly from the other nonperturbative effects. 

In contrast to the masses of u, d, and s quarks, masses of the heavy quarks must be 
retained in the hard-scattering terms in some cases. Several viable options exist for the 
treatment of {mg}, depending on the energy and fiavor composition of the scattering reac- 
tion. In two common approaches, the heavy quark contributions to the PDF's fa/A{x,fiF) 
and other resummed functions are either neglected or, alternatively, treated on the same 
footing as contributions from the light quarks (m, d, and s) above the respective heavy-quark 
mass thresholds, often placed, for convenience, at fip = rriQ. These choices correspond 
to the FFN and ZM-VFN factorization schemes, respectively. The FFN scheme optimally 
organizes the perturbative QCD series at energies of order of the heavy-quark mass ttiq. 
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The ZM-VFN scheme is the best choice for inclusive (depending on one momentum scale) 
distributions at energies much larger than tuq. It commonly utilizes dimensional regular- 
ization to expose collinear singularities of massless partonic matrix elements as poles 
in n = 4 — 2e dimensions. 

Neither of the two schemes is entirely satisfactory when applied to differential distribu- 
tions depending on two or more momentum scales of distinct magnitudes, as is the case of 
small-gT resummation at large Q. At sufficiently large \fS, the heavy-flavor quarks are copi- 
ously produced by quasi-collinear splittings of gluons along the directions of the initial-state 
hadrons. Such collinear contributions must be resummed at Q ^ uiq in the parton density 
fQ/A^x, fip) for the heavy quarks Q, so that a VFN factorization scheme is needed. If all 
scales (including qx) are of order Q, the heavy-quark mass rriQ <^ qx ^ Q can be neglected 
in the hard-scattering matrix elements, reducing the result to the ZM-VFN scheme. But 
niQ cannot be omitted at small qx {rriQ ^ qr Q), where it is not small compared to the 
momentum of the soft and collinear radiation. 

In 6-space, the form factor W{b,Q,XA,XB, {mQ}) is not well-deflned in the ZM-VFN 
scheme in the heavy-quark channels at impact parameters b > bo/rriQ, corresponding to 
factorization scales fip below the heavy-quark mass threshold, = b^/b < niQ. For charm 
quarks with mass nic = 1.3 GeV and bottom quarks with mass mt, = 4.5 GeV, the ZM-VFN 
form factor W is not well-defined at 6 > 0.86 and 0.25 GeV~\ respectively. The problem 
lies with the heavy-quark Wilson coefficient function CQ/a{x, b, uiq, fip), which is derived by 
using the factorization relation of Eq. (jH) with the densities Vq/a'{x, b, mo) and fa/A'{x, fip) 
of heavy quarks in partonic initial states A' = q, Q,g. At fip below the threshold, the fcr- 
integrated heavy-quark density fQ/A'{x, jJ^p) is set identically equal to zero, in accordance 
with the definition of fQ/A'{x, ^p) in the ZM-VFN scheme. Consequently the collinear poles 
1/e arising in the calculation of the 6-dependent density VQ/A'ixjb) for massless splittings 
Q ^ A' in 72 7^ 4 dimensions cannot be canceled by fa/A'{x, fJ'p) at such fip. Eq. ^ then 
implies that the Wilson coefficient function CQ/a{x,b,mQ, ftp) may be also infinite (contain 
the 1/e poles) below the heavy-quark threshold. The infinity arises because the ZM-VFN 
scheme incorrectly neglects the heavy-quark mass ttiq at energies fip = b^/b of order or less 
than mg. The solution to this problem is to retain the dependence on niQ at 6 > bo/niQ, 
which can be realized by formulating the CSS resummation in the GM-VFN factorization 
scheme. 
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The GM-VFN scheme [22| consistently implements the heavy-quark masses at all energy 
scales /iF, and it reproduces the FFN and ZM-VFN schemes in the limits ^ '^q and 
Uf ^ mo, respectively. Several versions of the GM-VFN scheme have been developed in the 
past years 0, H HQ Q Q, H H Q • 

A general procedure for the implementation of 
the GM-VFN scheme in qt resummation was outlined in Ref. where an application of 
the CSS resummation in the simplified Aivazis-Collins-Olness-Tung (S-ACOT) factorization 

nn 

scheme [2^ was presented for production of heavy quarks at HERA. 

The S-ACOT scheme is a variant of the GM-VFN scheme, which simplifies computations 
by neglecting the heavy-quark masses in hard subgraphs with incoming heavy quarks (flavor- 
excitation graphs), such as Q + g Q + 9- The heavy-quark masses are retained in the 
hard subgraphs with explicit production of heavy quarks (flavor-creation graphs), such as 
g + g Q+ Q. As demonstrated in Ref. ^ , application of the S-ACOT rules to the CSS 
resummed cross section efficiently retains the dependence on the heavy-quark masses where 
it is important, and drops it where it is not essential. 

In production of heavy gauge bosons {Q ^ w^q), the S-ACOT rules allow us to drop 
the rriQ dependence at qr ^ rriQ and keep the essential uiq dependence at qt ^ mQ. We 
assume that the heavy quarks are pairwise produced in perturbative splittings of gluons. 
We neglect possible, but yet experimentally unconfirmed, nonperturbative "intrinsic" heavy- 
quark contributions to the proton wavefunction Q]. We neglect the mass dependence 
entirely in the F-term, as it is non-negligible only at Qt ^ ttlq. In the M^-term, we drop 
the niQ dependence in all fiavor-excitation hard subgraphs and keep it in all fiavor-creation 
hard subgraphs. By this rule, rriQ is dropped in the perturbative Sudakov form factor S and 
the coefficient functions Ca/Q with the incoming heavy quarks, both of which are described 
by the fiavor-excitation Feynman graphs. We evaluate S up to O^a'^/n^) and Ca/Q up to 
0{as/iT) by using their massless expressions. 

We keep the ttiq dependence in the gluon-initiated coefficient functions Cq/q, since those 
are computed from the fiavor-creation Feynman graphs. The mass-dependent C(as/7r) co- 
efficient Cgj^ in Cg/g is given by 

7(1) 



CQ/g{x,b,mQ, ij,f) = Trx{1 - x)bmQKi{bmQ) 



KoibrriQ) - 6{fiF - mg) ln(^-^ ) 



(6) 



where -Pg/](a;) = Tr [x^ + (1 — xY] is the q ^ g splitting function, Tr = 1/2, and Kq{z) and 
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Figure 1: The 6-dependent parton densities Vqia{x, b, mg) vs. the impact parameter h for (a) charm 
quarks and (b) bottom quarks. The sohd and dashed curves correspond to the S-ACOT and massless 
("ZM-VFN") factorization schemes, respectively. 

Ki{z) are the modified Bessel functions [4^. The term proportional to the step function 
0{^p — mg) is nonzero above the heavy-quark threshold {^p > "^s), where the collinear 
logarithm P^^^^ (a;) In (/ii^/mg) is subtracted from CQjg{x,b,mQ, fif) and absorbed into the 
heavy-quark PDF fQ/A{x, fip)- The subtraction is not performed below the heavy-quark 
threshold, because fQ/A^x, fip) = at fip < uiq. For 6mg ^ 1, the gluon-initiated coeffi- 
cient function reduces to its massless expression, 



lim C^Ja;, 6, mg, /i^) = Trx{1 - x) - P^'/Jx) In {fipb/bo 



(1)/ 



b m 



(7) 



Convolutions of the Wilson coefficient functions CQ/a{x,b,mQ, fxp) and /cT-integrated 
PDF's fa/A{x,fiF) are combined into the 6-dependent parton density Pg/yi(x, 6, mg) ac- 
cording to Eq. (IH). At an n-th order of a^, Vq/a^x, b, mg) and its derivatives are continuous 
in ln(6) up to order n, if computed in the S-ACOT scheme. Here the order of PQ/^(a;, b, mg) 
is determined by a formal counting of powers of appearing both in the coefficient func- 
tions and PDF's. In this counting, the PDF fQ/A'{x, fip) for perturbatively generated heavy 
quarks is one order higher in than the light-parton PDF's. 

At leading order (LO), the function VQ/A{x,b,mQ) is composed of the heavy-quark- 
initiated coefficient function CQ/Q{x,b,mQ, fip), evaluated at order a^, and the gluon coef- 
ficient function CQ/g{x,b,mQ, fip), evaluated at order as/vr: 



VQ/Aix,b,mc 



C 



(0) 
Q/Q 



/qm) ix,b,mQ,HF) 
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+ (eg, ® /,m) (x,6,m2,/ii.) + 0(a^/vr^). (8) 

Here C^QiQ{x,h^mQ, hf) = 5(1 — x). The dependence of the leading-order (LO) parton 
density V Q/A{,x,h,mQ) on h is shown for the charm quarks in Fig. [ija) and bottom quarks 
in Fig. [ijb). The factorization scale /Xi? is set equal to h^/h to correctly resum the collinear 
logarithms in fQ/A{x, fip) in the small-6 limit (cf. Eq. Together with the S-ACOT 

predictions (solid lines), we show the ZM-VFN predictions (dashed lines). 

As discussed above, the ZM-VFN parton density V Q/Aix^h^mo) is not properly defined 
below the threshold iip = (or above b = bo/niQ). We formally define the "ZM-VFN" den- 
sity VQ/AiXjbjTriQ) at 6 > bo/niQ according to Eq. (jH)) by using massless Wilson coefficient 
functions Cj/a{x,b,mQ = Ojfip), as it was done in the previous resummation calculations. 
Such a definition provides just one of many possible continuations of the "ZM-VFN" coeffi- 
cient functions below the heavy-quark threshold, which render different results in qx space. 
We indicate this ambiguity by enclosing "ZM-VFN" in quotes and placing a question mark 
by the dashed line in Fig. [H At 6 > bo/niQ, the first term on the right-hand side of Eq. ((HI) 
vanishes (fQ/A{x, fip) = 0), so that VQ/A{x,b,mQ) becomes equal to the second term: 

"ZM-VFN", b > bo/niQ : 

Vq/a{x, b, ms = 0) = (eg, ® fg/A) {x, 6, ruQ = 0, fip) + Oia'Jn^). (9) 

The massless Q g coefficient function Cq/j,(x, 6, niQ = 0, fip) in Eq. ^ is given by Eq. 0. 

As in the ZM-VFN scheme, the heavy-quark PDF fg/A^x, fip) vanishes below the heavy- 
quark threshold in the S-ACOT scheme. However, the S-ACOT scheme properly preserves 
the mass-dependent terms below the threshold as a part of the gluon-initiated coefficient 
function Cgj,(x, 6, mg, /iir). The S-ACOT parton density 'Pq/a(x, 6, mg) is well-defined at 
all b. It reduces to the ZM-VFN result at 6 ^ bo/rriQ and is strongly suppressed at 6 3> 
bo/niQ by the modified Bessel functions KQ^bmo) and Kiipmo). The large-6 suppression 
is particularly strong in the case of the bottom quarks, where Vb/p{x,b,mi,) is essentially 
negligible at 6 > 1 GeV~^ (cf. Fig. [TJd). The suppression is caused by the decoupling of 
the heavy quarks in the parton densities at fip much smaller than ttiq {b much larger than 
bo /mo). This suppression is independent from the suppression of contributions with b > 
0.5 — 1 GeV~^ by the Sudakov exponential e~^^'''^^ (see Section ITl Djl . and it does not depend 
on Q. Consequently the impact of the nonperturbative contributions from b > 1 GeV~^ is 
reduced in the heavy-quark channels comparatively to the light-quark channels. 
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C. Nonperturbative contributions 



To gauge the effect of the nonperturbative contributions on the form factor W, we must 
estimate the behavior of the Sudakov exponential e"'^'^^''^) and ^-dependent parton densities 
Vj/A{x,b, {mQ}) at large b. Many studies have inve stigate d the nonperturbative contribu- 
tions to the resummed form factor; see, e.g., Refs. 0,l2a 22, 11, ^ M 



Here we opt to use a new parameterization of the nonperturbative component found in a 
global analysis Q| of Drell-Yan pair and Z boson production. This parameterization is 
obtained in the revised 6* model {3, 1^1, with the free model parameters chosen as to 
maximally preserve the exact form of the perturbative contributions at 6 < 1 GeV~^ 0|. 



The 6* model introduces a function b^{b,bmax) = b/ >/T+lPjl^^ and defines the re- 
summed form factor W as 

W{b, Q, XA, XB, {mQ}) = Wpert {b.{b, bma.), Q, xa, xb, {mg}) e'-^^^^^''^) (10) 

at all 6, with Wpert {b^{b,bmax),Q,XA,XB,{mQ}) being the finite-order ("perturbative") ap- 
proximation to W{b,Q,XA,XB, {itiq}). The higher-order corrections in and "power- 
suppressed" terms proportional to positive powers of b are cumulatively described by the 
nonperturbative function J^Np{b,Q), defined as 

^^P(6, Q) = - In f ^ W{b,Q,XA,XB,{mQ}) 

\Wpert {b*{b,b 

max ),Q,xa,xb, {mo}) 
and found by fitting to the data. 

In accordance with Ref. we choose bmax = 1-5 GeV~\ We also choose the factoriza- 
tion scale fip in (C /) {x, b, {ttiq}, fip) such that it stays above the initial scale Qo ~ 1 GeV 
for the PDF's fa/A{x,fiF)- Specifically, we define fip = bo/b^.{b,bo/Qo) = {Ql + bl/b'^y/'^, 
so that jjLp IS approximately equal to 6o/^ at 6 ^ and asymptotically approaches Qo at 
b ^ oo. 

The function J^Np{b, Q) found in the fit of Ref. 0| has the form 

J^Npib, Q) = 6^ [0.201 + 0.184 In (g/(3.2 GeV)) - 0.026 In (100 xaXb)] ■ (H) 

This parameterization has several advantages compared to the previous nonperturbative 
models. First, it minimizes modifications in Wpert{b,Q,XA,XB, {"^q}) by the 6* prescription 
in the small-6 region, where perturbation theory is valid. Second, Eq. (fTT| is in a good 
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agreement with both the available qt data and theoretical expectations. TNpip.Q) has a 
quadratic form, J-'np oc 6^, and follows a nearly linear dependence on ln(Q), which can be 
seen in Eq. (fTTll if the small term —0.026 ln(100 xa^b) is neglected. The quadratic form 
leads to Gaussian smearing by nonperturbative "intrinsic-ZcT" effects: W{h,Q) ~ e-o.5(fc|,>fe2^ 
with < 1 GeV^ in the observable Q range. Both the quadratic form and linear \n.{Q) 
dependence of J-'Npip,Q) are suggested by the infrared renormalon analysis of the large- 
h contributions The coefficient a2 = 0.184 GeV^ of the ln((5) term (found in the 

fit) agrees well with its independent estimate O.IQIqqq GeV^ obtained within the renor- 
malon analysis j^. The ln(Q) term arises from the Sudakov factor as a consequence of the 
renormalization-group invariance of the resummed form factor j^. It does not depend on 
the quark fiavors, and it contributes about 70% to the magnitude of J^Npip, Q) at Q Mz- 
Therefore, we expect the large fiavor-independent term 0.184 6^ ln(Q/3.2 GeV) to be also 
present in the heavy-quark scattering channels. 

The other terms in Eq. (fTT| may in principle depend on the fiavor of the participating 
quarks. However, no obvious dependence on the quark fiavor or nucleon isospin was ob- 
served in the global qt fit to the Drell-Yan pair and Z boson production data perhaps 
because these data are mostly sensitive to the scattering of light u and d quarks. Additional 
nonperturbative contributions may be present in the heavy-quark sector, especially if the 
"intrinsic" heavy-quark states 0| constitute a non-negligible part of the proton's wavefunc- 
tion. In the classical realization, the intrinsic heavy quarks contribute at large momentum 
fractions (x ^ 1), while the bulk of W and Z boson production happens at much smaller x, 
X = 10^^ — 10^^. Since the existence of the "intrinsic" heavy quarks remains an open ques- 
tion, we do not consider them in this study. We therefore parametrize the nonperturbative 
contributions by Eq. (fTI| . which neglects the fiavor dependence. We vary the parameters of 
^Np{b,Q) in the heavy-quark channels in order to evaluate the uncertainty in the physical 
observables arising from this approximation. 



D. An example of the 6-space form factor 

Fig. [21 shows an example of the full resummed form factor bW{b,Q) (cf. Eq. (jSJ) 
in Z boson production via bb annihilation in the Tevatron Run-2. The form factor is 
computed in the S-ACOT scheme at the leading-order (LO) and in the S-ACOT and 
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p p ^ Z° X (VS = 1 .96 TeV); Q=91 .1 87 GeV; y=0; CTEQ5HQ1 
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Figure 2: The form factor bW{b,Q) vs. the impact parameter b for bb — > Z^X in the Tevatron 
Run-2. The sohd, dashed, an dotted curves correspond to the NLO S-ACOT, NLO "ZM-VFN", and 
LO S-ACOT calculations, respectively. 

"ZM-VFN" schemes at the partial next-to-leading-order (NLO). The leading-order "ZM- 
VFN" contribution is computed by using the formal definition of Pq/^(x, 6, mg = 0) 
in Eq. (jH]). The NLO curves are obtained by adding the massless flavor-excitation con- 
tribution (^g^/Q® Iq/aJ {x,b,mQ = 0,/iir) of order C(as/7r) to the LO parton density 
Pq/a(x, 6, mg) in Eq. We do not include the O^al/n^) flavor-creation contribution 

i^Q/a ® fa/A^ {x,b,mQ, fip), vi^ith a = q, g, needed to evaluate VQ/A{x,b, {niQ}) at the full 
NLO accuracy. The partial NLO approximation results in a discontinuity of the derivative 
dW{b,Q)/d\nb at the bottom quark threshold, which has little impact on our results. The 
perturbative contribution to the Sudakov factor S{b,Q) is evaluated up to order a^/vr^. 
The nonperturbative terms are included according to Eqs. (fTn|l and (fTI| . assuming flavor 
independence of J-'Npib,Q). Furthermore, the CTEQ5HQ1 parton distributions |^ were 
used. 

According to the figure, the form factor bW(b,Q) in bb ^ Z exhibits a maximum at 
b ~ 0.07 GeV~\ Larger impact parameters are suppressed by the Sudakov exponential 
g-5(6,Q) amount growing with Q, independently of the quark flavor and factorization 

scheme. In the S-ACOT scheme, the form factor is also suppressed by the mg decoupling 
in V Q/A{,x,b,mQ). Consequently the mass-dependent variations in bW{b,Q) are more pro- 
nounced in production of not too heavy bosons, such as Z^, and in bottom quark channels, 
where they occur at relatively small b ~ 0.25 GeV~^ (cf. Fig. [l]). 
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The "ZM-VFN" form factor underestimates the S-ACOT form factor in a wide range 
of b, and it exhibits non-smooth behavior in the vicinity of the bottom quark threshold. 
The NLO S-ACOT and "ZM-VFN" form factors agree well at 6 < 0.05 GeV"\ The NLO 
correction to Vq/a^x, b, niQ) mainly contributes at 6 ~ 0.02 — 0.25 GeV~\ where it enhances 
the LO form factor. The NLO correction is turned off at 6 > 0.25 GeV~^ by the condition 
/q/^(x, /ii?) = for hf < fTiQ. It is also suppressed at 6 < 0.02 GeV~^ by the small as{bo/b). 

The differences between the S-ACOT and "ZM-VFN" form factors seen at 6 ^ bo/mQ will 
affect the cross sections in qr space at small and moderate transverse momenta. 



III. NUMERICAL RESULTS 



In the present section we compare the resummed qt distributions calculated in the massive 
(S-ACOT) and massless ("ZM-VFN") schemes, defined according to Eq. We examine 
Drell-Yan production of W^, and Higgs bosons in the Tevatron Run-2 (the cm. energy 
= 1.96 TeV) and at the LHC {y/S = 14 TeV). In the case of W boson production, we 
consider the leptonic decay mode W ^ eu and discuss the impact of the different schemes on 
the measurement of the W boson mass. We apply the resummation formalism described in 
the previous section. The numerical calculation was performed using the programs Legacy 
and ResBos |43, with the CTEQ5HQ1 parton distribution functions We use the 
parameterization (|TT|l for the nonperturbative function J-'Np{b,Q), unless stated otherwise. 
The and Z° boson masses are assumed to be Mw = 80.423 GeV and Mz = 91.187 GeV, 
respectively. The heavy quark masses are taken to be rric = 1.3 GeV and nih = 4.5 GeV. 



A. Partonic subprocesses in W and Z" boson production 

We first compare contributions of various partonic subprocesses to the total and Z° 
production cross sections. These contributions will be classified according to the types of 
the quarks q and antiquarks q' entering the qq'W or qqZ electroweak vertex. At the Born 
level, the cross section for production of narrow-width and Z° bosons in a qq' partonic 
channel is approximately given by 

(^qq' ~ CTQ ^qq' (t) dgq' , 
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W+ 


Subprocesses 


ud 


us 


cd 


cs 


c6 


Tevatron Run-2 


90 


2 


1 


7 





LHC 


74 


4 


1 


21 






w- 


du 




dc 


sc 


6c 


90 


2 


1 


7 





67 


2 


3 


28 








uu 




ss 


cc 


66 


57 


35 


5 


2 


1 


36 


34 


15 


9 


6 



Table I: Partial contributions a^-,/atot of quark-antiquark annihilation subprocesses to the total 
Born cross sections in and boson production at the Tevatron and LHC (in percent). 



where r = Q'^/S, Q = Mw {Mz) in (Z°) boson production, 



/,/A(e, Q) U/b{j, Q) + /,-'M(e, Q) U/b{j. Q) 



(12) 



is the parton luminosity in the channel, is the flavor-dependent prefactor, and a. 
includes the remaining flavor-independent terms |2J|. In boson production, 

2 



9, 



qq' 



V- 
V' 



where Vg-i is the appropriate CKM matrix entry. In boson production, 

9qq' = 5qq' [(1 - 4 Slu^ O^)"^ + l] , 



(13) 



(14) 



where = 2/3 or —1/3 is the quark electric charge in units of the positron charge, and 
is the weak mixing angle. The partial cross sections aqq,/(Jtot (where atot = J2qq'^qq') 
in the different qq' channels are evaluated as C-qqig^qi/ {^qqi '^qq'Qqq'^ and listed in Table HI 
These values will be modified somewhat by the NLO radiative corrections included in the 
following subsections. 

According to the Table, the contributions of the heavy-quark channels at the Tevatron are 
small. In W boson production, only 8% of the cross section involves heavy (charm) quarks, 
while light quarks account for 92%. In Z boson production, heavy quarks contribute even 
less, with only 3%. Due to the small fractional contribution of the heavy quarks at the 
Tevatron, the differences between the massive and the massless calculations are negligible 
given the expected experimental uncertainties. 

The LHC probes smaller momentum values x, where the fractional contributions of the 
initial charm and bottom quarks are much larger. These channels sum up to 22% in bo- 
son production, 31% in W~ boson production, and 15% in Z^ boson production. Therefore, 
the following discussion will concentrate on the LHC and only briefiy mention the Tevatron 
processes. 
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[GeV] [GeV] 

(a) (b) 
Figure 3: da/dgx for cc, bb boson production at the LHC: (a) bb channel only, (b) combined 

cc and bb channels. The solid (red) curve shows the distribution in the massive (S-ACOT) scheme. 
The dashed (black) curve shows the distribution in the massless ("ZM-VFN") scheme, computed 
using the parametrization ifTTIl of the nonperturbative function J^Np{b,Q). The dot-dashed (blue) 
line was calculated in the "ZM-VFN" scheme using an alternative parameterization Q| of the 
nonperturbative function J^Np{b,Q)- 

B. boson production at LHC 

We first examine the infiuence of the charm and bottom quark masses on the resummed 
Z° boson cross section at the LHC. Fig. El compares the transverse momentum distribution 
of boson production in the heavy-quark (cc and bb) channels at the LHC, calculated in 
the S-ACOT and "ZM-VFN" schemes. Fig. [H^a) shows the distribution in the bb channel 
only, where the effect is the largest due to the large b quark mass. In the approximate "ZM- 
VFN" scheme (black dashed line), the da/dgx distribution is shifted in the maximum region 
at ~ 13 GeV to larger transverse momenta by about 2 GeV. This reflects the enhancement 
of the form factor b W{b, Q) at b ~ bo/nib ^ 0.25 GeV"^ in the S-ACOT scheme. The region 
of large transverse momenta (g^ ^ 25 GeV) is essentially sensitive only to small impact 
parameters {b < 0.1 GeV^^), where the S-ACOT and ZM-VFN form factors are very close. 
Consequently, the two schemes give very close predictions in the large-gr region. 
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The shift to the larger values of qr is similar in the cc channel, albeit smaller due to the 
lower charm mass. Fig.E^b) shows the combined distribution in the cc and bb channels. The 
mass-induced shift of the peak of the qt distribution is much smaller, but clearly visible. One 
should notice that the distribution of Fig.lHl^b) peaks around 7.5 GeV, since it is dominated 
by the charm contribution, whereas the bottom distribution peaks at 15 GeV (Fig. [Hfa)). 
The combined effect, displayed in Fig.lHl^b), is an enhancement of the rate in the peak region 
of about 2.5%. 

For comparison, the dot-dashed (blue) line shows the shift of the "ZM-VFN" qr dis- 
tribution due to the choice of a different parameterization for the nonperturbative func- 
tion J^Np{b, Q) j^l- The difference between the dotted (black) line and the dot-dashed (blue) 
line provides a conservative estimate of the current experimental uncertainty in J^Np{b,Q). 
The mass-induced effect is at least comparable to the other uncertainties, when the cc and bb 
channels are taken alone. However, these channels cannot be isolated from the light-quark 
channels in the single-particle inclusive observables. According to Table HI the heavy-quark 
channels contribute about 15% to the total Z° boson cross section. If the light-quark flavors 
are included, the S-ACOT rate is enhanced in the maximum region by a small value of 
~ 0.3%, which is beyond the experimental precision of the LHC. 



C. W boson production at LHC 

boson production at the LHC is of particular interest, as it will be utilized to measure 
the W boson mass M^y with high precision (tentatively less than 15 MeV). The W boson 
mass can be extracted in the leptonic decay W ^ eu hj fitting the theoretical model for 
different values of Mw to kinematical Jacobian peaks, arising at M^" ^ My/ in the leptonic 
transverse mass [M^] distribution and at ~ Miy/2 ^ 40 GeV in the electron 
spectrum. 

The S-ACOT and "ZM-VFN" schemes yield essentially identical results for the M^"" dis- 
tribution, suggesting that the error caused by the massless approximation is likely negligible 
in the method. This conclusion follows from the low sensitivity of the distribution 
to variations 5qT in the bosonic qx distributions (suppressed by 6q^/M^ <^ 1). However, 
the distribution may be strongly affected. 

To better display percent-level changes in da/dp^ associated with the effects of the heavy 
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(a) (b) 
Figure 4: The fractional change in the distribution da/dp^ for pp W^X e^V(.X at the 
LHC: a) only partonic channels containing heavy quarks {ub,cd,cs,cb W~) are included; (b) all 
channels are included. The dashed (black) line shows the relative difference between the massless 
("ZM-VFN") and S-ACOT cross sections. The dotted (green) and dot-dashed (blue) lines show the 
relative changes in the S-ACOT cross section induced by variations of M\y by ±30 MeV. 

quarks, Fig. IH shows the fractional difference [da"'°'^ / dp^) / {da^^'^^^'^ / dp'j.^ — 1 of the 
"modified" (e.g. "ZM-VFN") cross section and the "standard" (S-ACOT) cross section. We 
compare the modifications due to the approximate massless "ZM-VFN" treatment to the 
modifications due to explicit variations of in the S-ACOT result. 

We first examine the partonic production channels with at least one c or 6 quark (or their 
antiparticles) in the initial state. Fig. Ufa) shows the fractional change of the p^ distribution 
for these channels in W~ boson production. The solid red line shows the reference result, 
calculated using the S-ACOT scheme with Mw = 80.423 GeV. The dotted (green) and dot- 
dashed (blue) lines show the variation of the S-ACOT p^ distribution, if the W boson mass 
is shifted by ±30 MeV. A small increase of Mw results in a positive shift of the Jacobian 
peak (dotted green line), which reduces the cross section at pf. < Mw/'^ and increases the 
cross section at > M^jl. Smaller values of M^J^ result in a shift in the opposite direction. 

The dashed black curve is the ratio of the p\. distributions calculated in the "ZM-VFN" 
and S-ACOT schemes for = 80.423 GeV. In the heavy-quark channels only (Fig. (jlfa)). 
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Figure 5: Rapidity distributions of W bosons at the LHC. The solid (dot-dashed) line is the 
distribution of the bosons in all channels. The dashed (dotted) line is the distribution of 

the Ty+(iy~) bosons in the channels that involve at least one initial-state c or 6 quark. 

the "ZM-VFN" approximation shifts the Jacobian peak in the positive direction. The size 
of this effect is comparable to a shift in Mw of about +35 MeV. 

Fig.EJa) is also representative of boson production at the LHC, because the dominant 
heavy-quark contributions (cs and cs — > W~) are charge-conjugated, and the parton 

density functions of c (s)-quarks and c (s) -anti-quarks are almost identical. Hence, the rates 
for and W~ boson production in heavy-quark scattering are very similar. 

Fig. mb) shows the fractional change in the distribution for W~ bosons, summed 
over all possible partonic states. The total shift of the differential cross section due to the 
"ZM-VFN" approximation is comparable with a positive My/ shift of about 10 MeV. This 
value is consistent with the shift of 35 MeV in the c and h channels, in view of the fact that 
the c, 6-channels contribute ~ 30% to the total W~ cross section according to Table HI The 
qT shift is independent of the rapidity y of the W~ boson, because the rapidity distributions 
of W~ bosons have similar shapes in the heavy-quark and combined channels (see Fig. [HI). 
Consequently the mass-induced shift in the W~ boson distribution is not affected by the 
rapidity cuts. 

We contrast this result with boson production, where the profile of the rapidity 
distribution for the heavy-quark channels differs from that for the combined channels. The 
rapidity distribution of bosons in all channels (solid line in Fig. EJ has characteristic 
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shoulders at ?/ ~ ±3. These shoulders are produced by ud,us scattering, enhanced 

at large momentum fractions x by the valence w-quark distribution. The fraction of c and 
h channels (dashed line) is the largest at central rapidities, — 1 < y < 1. In this region, c 
and h contribute ~ 29% to the cross section. The difference between the "ZM-VFN" and 
S-ACOT distributions is comparable here to a change in Mw of about 9 MeV. In the 
full rapidity range, this difference is comparable to 5Mw ^ 6 MeV. If experimental cuts 
appropriate for the ATLAS experiment (|?/e| < 2.5, > 25 GeV and > 25 GeV) are 
applied, the difference is comparable to 6Mw ^ 7.5 MeV. 

The presented estimates may be modified by additional nonperturbative contributions 
arising in the heavy-quark channels, such as those associated with the "intrinsic" heavy 
quarks. Our model has neglected such contributions, in view that their magnitude is un- 
certain, and they are less likely to contribute at the small x typical for W and Z boson 
production. We assumed that the dominant nonperturbative contribution at Q ~ Mz 
(~ 70% of J-'Np{b,Q)) arises from Q dependence of the Sudakov factor, which does not 
depend on the quark fiavor (cf. Section Hi Cjl . The "ZM-VFN" and S-ACOT functions differ 
in b space by several tens of percent in the threshold region (cf. Figs. [U and EJ. The ex- 
tra nonperturbative terms must contribute at a comparable level to be non-negligible. We 
confirmed this estimate numerically by repeating the analysis of the distributions for a 
varied (rescaled by a factor of two) function J^Npip.Q) in the heavy-quark channels. The 
resulting variation in pf. distribution for the heavy-quark channels, comparable with 5Mw 
of a few tens MeV (???), is consistent with the effect of the rescaling of J-'Npip, Q) in h space, 
of order of the difference of the "ZM-VFN" and S-ACOT form factors W{h, Q) in the large-6 
region. 

D. Z° and boson production at the Tevatron 

The dependence on the heavy-quark masses in and boson production at the 
Tevatron follows nearly the same pattern as that at the LHC. In Z^ boson production via c 
and h quark channels at the Tevatron, the qx distribution in the "ZM-VFN" approximation is 
shifted toward larger g-r with respect to the S-ACOT distribution, in analogy to the similar 
shift at the LHC (cf. Fig.Oj). The maximum of the qr distribution in the hh channel is shifted 
by approximately 2.5 GeV if the masses of the quarks are neglected; this is slightly larger 
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than the shift at the LHC. However, modifications in the qt distribution in the combined 
initial-state channels are negligible, because c and b scattering contributes only 3% to the 
total cross section (cf. Table P). 

The analysis of W boson production at the Tevatron is less involved, given that the rates 
for the processes pp — ^ W~^X and pp — ^ W~X are related by charge conjugation. The 
electron p^ distributions in the c and b quark channels at the Tevatron qualitatively follow 
the pattern shown in Fig. HJa) for the LHC. The difference between the p^ distributions 
in the "ZM-VFN" and S-ACOT schemes at the Tevatron is comparable to the effect of a 
Mw increase by 45 MeV (vs. 35 MeV at the LHC). The heavy-quark channels contribute 
just 8% to the total cross section at the Tevatron. Consequently the net Mw shift in 
all channels (3 MeV in the p^ method and in the M^" method) is tiny compared to the 
expected experimental uncertainty of ~ 30 — 40 MeV. 



E. Higgs boson production 

In light of our analysis of the bb Z° channel fSec. lIIIB|l . we anticipate that production 
of Higgs bosons via bb annihilation is also sensitive to the heavy-quark mass effects. 
In the standard model, the gluon-gluon fusion process gg is the dominant produc- 

tion channel for a wide range of Higgs boson masses at both the Tevatron and LHC (see, 
e.g., 0,y, Ol and references therein). Following the gluon-gluon fusion, the next most impor- 
tant channels at the Tevatron are associated WH^ and ZH^ production, and vector boson 
fusion. At the LHC, the next largest production modes are the electroweak boson fusion 
and Higgs-strahlung processes. The process bb — > typically contributes less than 1% of 
the total production cross section. In view of the current uncertainty in the Higgs transverse 
momentum distribution effects due to the heavy quark mass in the bb channel will be 
negligible in the standard model. 

The question becomes more interesting in extensions of the standard model containing 
two Higgs doublets, like the minimal supersymmetric standard model (MSSM), where the 
Yukawa couplings of the bottom quarks to the neutral Higgs bosons h^, and ^4° depend 
on the supersymmetric parameter tan /3. If tan /3 is large, the bottom Yukawa coupling 
is strongly enhanced, while the top quark Yukawa coupling remains nearly constant or is 
suppressed (see Ref. 0|). The bb annihilation rate is comparable to the gluon-gluon fusion 
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Figure 6: Transverse momentum distribution of on-shell Higgs bosons in the hh — > channel 
at (a) the Tevatron and (b) LHC. The soHd (red) Hues show the qT distribution in the massive 
(S-ACOT) scheme. The dashed (black) hnes show the distribution in the massless ("ZM-VFN") 
scheme. 

rate for medium values of tan/? and can dominate the cross section for tan/5 > 30. 

We implemented the heavy-quark mass effects in the resummation subroutine for 66 — > if ° 
developed in an earlier study |^.^ Fig. IHl displays da/dqT for hb — > boson production at 
(a) the Tevatron and (b) LHC. The Higgs boson mass is chosen to be Mh = 120 GeV. The 
supersymmetric bb — Higgs couplings are obtained at leading order bvrescaling the standard 
model coupling: 9bi{h^'^H0 a^} ~ {— sin a, cosa, sin/375}(7^^/cos/? 0, 0|. The net effect 
of the bottom quark masses on qt distributions will be the same for both SM and MSSM 
neutral Higgs bosons, up to an overall normalization constant. For this reason, Fig. IHl does 
not specify the overall normalization of Qt distributions. 

The S-ACOT and "ZM-VFN" distributions da/dgr are shown by the solid (red) and 
dashed (black) lines, respectively. A significant shift of the distribution to larger values 
of qt is seen in the "ZM-VFN" approximation. As in bb process (Sec. IIII B|) . the 

enhancement of the qt distribution at small qt is caused by the enhancement in the S- 
ACOT form factor bW{b, Q) at intermediate and large b {b> 0.1 GeV~^). At the Tevatron, 
Fig. Efa), the qr maximum shifts in the "ZM-VFN" approximation to larger by about 
2 GeV out of 11.7 GeV (about 17%). For a Higgs mass Mu = 200 GeV, the maximum of 



^ We thank A. Belyaev and C.-P. Yuan for pointing out a typo in the C^/q function in Ref. 0|, which does 
not affect the results shown in this section. 
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Differences in da / dqx for hb — > at the LHC 



mu [GeV] 


120 


250 


600 


Position of the 


"ZM-VFN" 


15.4 


16.8 


18.8 


maximum [GeV] 


S-ACOT 


14.1 


15.8 


18.2 


Difference in the positions [GeV] 


1.3 


1.0 


0.6 



Table II: Positions of the maximum of the distribution da/dqx in bb — > process at the LHC in 
the mass-dependent (S-ACOT) and massless ("ZM-VFN") calculations, as well as the differences of 
these positions. 

da/dgx shifts by about 1.9 GeV out of 12.7 GeV. 

The corresponding qt distributions for the LHC are shown in Fig. Efb). The difference 
between the "ZM-VFN" and S-ACOT calculations is smaller compared to the Tevatron, 
because the influence of the region b > 0.1 GeV~^ is reduced at smaller momentum fractions x 
probed at the LHC 0|. The maximum of the transverse momentum distribution shifts in 
the "ZM-VFN" approximation by about 1.3 GeV (9% out of 14.1 GeV) to larger qt- The 
results for two other Higgs masses {Mh = 250 and 600 GeV) are summarized in Table ITTl 
The full qr dependence of the bb H process is affected by additional factors, such as the 
constraints on phase space available for QCD radiation (less relevant at small gr)- The full 
qt dependence is investigated elsewhere |6l| . 



IV. CONCLUSION 



In the present paper, we make use of recent theoretical developments to estimate the 
impact of heavy-quark masses mg on transverse momentum (qt) distributions in production 
of W, Z, and Higgs bosons at the Tevatron and LHC. We note that the the zero-mass 
("massless") variable flavor number scheme is not consistent in the heavy-quark channels, 
if qt is of order of the heavy-quark mass uiq. To properly assess the mg dependence, 
we perform resummation of the large logarithms ln(g|n/Q^) in the Collins-Soper-Sterman 
approach, formulated in a general variable flavor number factorization scheme (S-ACOT 
scheme) to correctly resum the collinear logarithms ln(/Xir/mg). 



22 



We compare our consistent treatment of the heavy-quark mass dependence with an ap- 
proach based on the ZM-VFN scheme. The proper treatment of the heavy-quark masses 
leads to an enhancement of the form factor W{b, Q) at intermediate and large impact pa- 
rameters, which in turn increases the transverse momentum distribution at small qr- In 
Drell-Yan production of heavy bosons, the cumulative effect in the S-ACOT scheme shifts 
the peak of the da/dqr distribution to smaller values of qr. The mass effects are negligible 
if qr is of order Q. 

The mass dependence is magnified in subprocesses dominated by scattering of initial- 
state bottom quarks, as a consequence of the larger 6-quark mass. In the production of Z° 
and light Higgs bosons via hb annihilation at the LHC, the maximum of da/ dgr shifts from 
13.5 - 14 GeV in the S-ACOT scheme to ^ 15 GeV in the "ZM-VFN" approximation. At 
the Tevatron, the maximum shifts from ~ 10 to 12 GeV. Our computation is applied to 
single-particle inclusive production, when no heavy quarks are observed in the final state. 
Since the single-particle inclusive cross sections receive contributions from both the heavy- 
quark channels and intensive light-quark channels, the uiq dependence is reduced below the 
experimental sensitivity in several cases, notably in and boson production at the 
Tevatron, and Z° boson production at the LHC. On the other hand, the signal from the 
heavy-quark channels may be enhanced by tagging one or two heavy quarks in the final 
state. In that case, mass-dependent effects of a similar magnitude may be observed in some 
regions of phase space, e.g., at small qt of the boson-heavy quark system. 

Since the LHC plans to measure the W boson mass Mw to a high precision, we must fully 
account for all potential uncertainties, including the heavy-quark mass effects. At the LHC, 
the massless "ZM-VFN" approximation introduces a bias in the measurement of Mi^ from 
the transverse momentum distribution da/dp^ of the electrons, comparable to a positive 
MiY shift by ~ 10 MeV. In boson production the bias is comparable to the shift of 
9 MeV at central rapidities of W bosons (| y | < 1) and 6 MeV in the whole rapidity range. 
These effects are not negligible in view of the desired precision of the Mw measurement of 
~ 15 MeV. 
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